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The Boltzmann equation for inelastic Maxwell models is considered in order to investigate the 
dynamics of an impurity (or intruder) immersed in a granular gas driven by a uniform shear flow. 
The analysis is based on an exact solution of the Boltzmann equation for a granular binary mixture. 
It applies for conditions arbitrarily far from equilibrium (arbitrary values of the shear rate a) and 
for arbitrary values of the parameters of the mixture (particle masses rm, mole fractions Xi, and 
coefficients of restitution Qij). In the tracer limit where the mole fraction of the intruder species 
vanishes, a non equilibrium phase transition takes place. We thereby identity ordered phases where 
the intruder bears a finite contribution to the properties of the mixture, in a region of parameter 
space that is worked out in detail. These findings extend previous results obtained for ordinary 
Maxwell gases, and further show that dissipation leads to new ordered phases. 

PACS numbers: 05.20.Dd, 45.70.Mg, 51.10.+y 



I. INTRODUCTION 

Inelastic hard spheres (IHS) provide a useful theo- 
retical, computational, and experimental framework for 
studying granular gases [1-3]. In the simplest version, 
the spheres are assumed to be smooth (i.e. frictionless) 
and the inelasticity in collisions is specified in terms of a 
constant coefficient of normal restitution a < 1 [1] . At a 
kinetic theory level, the essential information about the 
dynamical properties of the fluid is embedded in the one- 
particle velocity distribution function. For sufficiently 
low-densities, the conventional Boltzmann equation can 
be extended to IHS by changing the collision rules to 
account for the inelastic character of the interactions 
[4, 5]. However, the complex mathematical structure 
of the Boltzmann collision operator prevents one from 
obtaining exact results, so that most of the analytical 
progress for IHS has been achieved by using approxi- 
mated methods and/or simple kinetic models [6]. For 
instance, the explicit forms of the Navier-Stokes trans- 
port coefficients have been approximately determined by 
considering the leading terms in a Sonine polynomial ex- 
pansion [7]. 

Needless to say, the difficulties of solving the Boltz- 
mann equation increase considerably when one considers 
far from equilibrium situations, such as shear flow prob- 
lems. Although good estimates for nonlinear transport 
properties of IHS have been obtained [8], the search of 
exact results has stimulated the use of model interactions 
simpler than hard spheres. One possibility is to consider 
a mean field version of the hard sphere system, where 
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randomly chosen pairs of particles inelastically interact 
with a random impact direction. This assumption yields 
a Boltzmann equation where the collision rate is indepen- 
dent of the relative velocity of the two colliding particles 
[9]. This interaction model is usually referred to as the 
inelastic Maxwell model (IMM). It must be noted that 
in the conventional case of ordinary gases colliding elasti- 
cally, Maxwell models correspond -in three dimcnsions- 
to particles interacting via a repulsive potential propor- 
tional to the inverse fourth power of distance [10]. Nev- 
ertheless, in the framework of the Boltzmann equation, 
one can introduce Maxwell models at the level of the 
cross section, without any reference to a specific inter- 
action potential [11, 12]. Thanks to the simplifications 
introduced by IMM in the kernel of the Boltzmann colli- 
sion operator, it is possible in some particular situations 
to find non-trivial exact solutions to the Boltzmann equa- 
tion [13-18]. Apart from their academic interest, it has 
also been shown in some cases that the results derived 
from IMM [13] agree well with those obtained analyti- 
cally for IHS from Grad's method [19, 20] and by means 
of Monte Garlo simulations [19]. We will substantiate 
this point in the concluding section. In addition, even 
recent experiments [21] for magnetic grains with dipolar 
interactions are qualitatively well described by IMM. All 
of these results clearly show the utility of IMM as a toy 
model to unveil in a clean way the influence of the in- 
elasticity of collisions in granular flows, especially in far 
from equilibrium situations where simple intuition is not 
enough. 

The aim of this paper is to investigate the dynam- 
ics of an intruder or impurity immersed in an inelastic 
Maxwell gas subject to the simple or uniform shear flow 
(USF). This state is perhaps one of the most widely stud- 
ied states in granular gases [22]. From a macroscopic 
point of view, the USF is characterized by constant par- 
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tial densities uniform granular temperature T, and 
a linear velocity field = ay, a being the constant shear 
rate. Because the only hydrodynamic gradient is that of 
flow velocity, the mass and heat fluxes vanish by sym- 
metry and the pressure or stress tensor Pij is the only 
relevant flux in the problem. The knowledge of the ele- 
ments of Pij gives access to the rhcological properties of 
the mixture. In the context of IMM, the above transport 
properties have been recently obtained [18] in terms of 
the shear rate and the parameters of the mixture. Our 
goal now is to consider the tracer limit. Since this case 
corresponds to a situation in which the mole fraction 
Xi = ni/(ni + n2) of the tracer species is negligible, one 
expects that the properties of the excess species (granu- 
lar gas) are not affected by the presence of the impurity 
particle. In particular, the relative contribution of the 
impurity to the total energy of the system, Ei/ E, is ex- 
pected to be proportional to xi when xi <^ 1 (here, E 
denotes the total kinetic energy of the system and Ei is 
the kinetic energy of the impurity). Consequently, the 
contribution of the impurity to the total energy is likely 
to be negligible. However, as in the elastic case [23], we 
present in this paper a violation of the above expectation, 
that is ascribable to a non-equilibrium phase transition. 
Specifically, we found that the impurity particle has a 
finite contribution to the total energy of the system (im- 
purity plus granular gas) either for light impurities under 
conditions that depend critically on the asymmetry of 
collisional dissipation (impurity-gas against gas-gas col- 
lisions), or conversely, for heavy impurities provided the 
shear rates are not too large. This phenomenon is akin 
to a phase transition, where the order parameter is given 
by the ratio of the impurity kinetic energy over the to- 
tal kinetic energy of the system. We borrow here the 
terminology of equilibrium phase transitions; it should 
nevertheless be stressed that the kind of ordering under 
discussion is not spatial, but refers to a specific, kinetic, 
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Here, 

= y rfv/r(v) (3) 

is the number density of species r, w^s is an effective 
collision frequency (to be chosen later) for collisions of 
type r-s, = 27r'^/^/r(d/2) is the total sohd angle in d 
dimensions, and a^-g < 1 refers to the constant coefficient 
of restitution for collisions between particles of species r 
with s. In addition, the primes on the velocities denote 
the initial values {v'j^,V2} that lead to {vi,V2} following 
a binary collision; 

v'l = vi - (l + a,T/) (ct • gi2)ff, (4) 



order parameter (Ei/E); as such, the "ordering" sce- 
nario specifically belongs in non-equilibrium. Moreover, 
we found that ordering always sets in at large enough 
shear rates (i.e., for shear rates larger than a certain crit- 
ical value), provided the impurity is sufficiently light com- 
pared to the grains of the host gas. While this feature is 
common to elastic [23] and inelastic systems, the possi- 
bility of what we refer below as heavy tracer ordering, is 
specific to inelastic gases. A preliminary account of part 
of this work has appeared in Ref . [24] . 

The plan of the paper is as follows. The Boltzmann 
equation for IMM is introduced in Sec. II and the col- 
lisional moments needed to get the pressure tensor are 
explicitly evaluated. In Sec. HI, the USF problem is 
presented and the rheological properties are obtained in 
terms of the shear rate, the masses, the concentration, 
and the coefficients of restitution. The main results are 
derived in Sec. IV. Specifically, the tracer limit is con- 
sidered in detail, which shows the existence of the above 
non-equilibrium transition. Finally, the paper is closed 
in Sec. V with a brief discussion. 

II. THE INELASTIC MAXWELL MODEL 

Let us consider a binary mixture of inelastic Maxwell 
gases at low density. In the absence of external forces, 
the set of nonlinear Boltzmann equations for the mixture 
reads 

(^+v-v)/.(r,v;t)=^J,,[vl/,(t),/,(t)] , (1) 



where fr{r,v;t) is the one-particle distribution function 
of species r (r = 1,2) and the Boltzmann collision op- 
erator Jrs [vl/r,/s] describing the scattering of pairs of 
particles is 



(2) 

I 

V2 = V2 + ^-lrs (1 + a^/) {a ■ gia)? , (5) 

where gi2 = vi— V2 is the relative velocity of the colliding 
pair, is a unit vector directed along the centers of the 
two colliding spheres, and firs = nir/{mr+ms). From the 
densities n^, we define the mole fractions Xr = nr/{ni + 
"-2)- 

The effective collision frequencies cj^s are independent 
of the relative velocities of the colliding particles but de- 
pend in general on space and time through its depen- 
dence on density and temperature. In previous works on 
multicomponent granular systems [13, 18, 25], ujrs was 
chosen to guarantee that the cooling rate for IMM be 
the same as that of the IHS. With this choice ("improved 



Jrs[vi\fr,fs] = :^^Jdv2 J da [a^} fr{r , v[ , t) fs{r , v'^, t) - fr{r , vi , t) fs{r , V2 , t)] . 
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Maxwell model"), the collision rates oJrs are functions of 
the temperature ratio 7 = T1/T2, that is itself a func- 
tion of the (reduced) shear rate in the USF problem. 
A consequence of this choice is that one has to numeri- 
cally solve a set of nonlinear equations in order to get the 
shear rate dependence of the temperature ratio [13, 18]. 
Since we wish to obtain in this paper analytical results 
for arbitrary spatial dimensions in a quite complex prob- 
lem that involves a delicate tracer limit, we will consider 
here a simple version of IMM where w^s is independent 
of the partial temperatures of each species ( "plain vanilla 
Maxwell model"). Thus, one defines UJr a as 



and the heat flux 



An, 



(6) 



where the value of the constant A is irrelevant for our 
purposes. Here, n = rir is the total number density 
of the mixture. The form of oJrs is closer to the original 
model of Maxwell molecules for elastic mixtures. This 
plain vanilla model has been previously used by several 
authors [26-29] in some problems pertaining to granular 
mixtures, and we will argue in section V that it is capable 
of capturing the essential physical effects at work here. 

At a liydrodynamic level, the relevant quantities in a 
binary mixture are, apart from n^, the flow velocity u, 
and the "granular" temperature T . They are defined in 
terms of moments of the distribution fr as 



pu : 



nT 



E 



rfv^FV.(v), 
a 



(7) 



(8) 



(11) 



respectively. Finally, the rate of energy dissipated due to 
collisions among all species defines the cooling rate C as 



J d\-V'^Jrs[v\fr, fs] = -dnT( . 



(12) 



The key advantage of the Boltzmann equation for 
Maxwell models (both elastic and inelastic) is that the 
(coUisional) moments of Jrsifn fs] can be exactly evalu- 
ated in terms of the moments of fr and fs without the 
explicit knowledge of both velocity distribution functions 
[34, 35] . This property has been recently exploited [36] to 
obtain the detailed expressions for all the second-, third- 
and fourth-degree coUisional moments for a monodisperse 
gas. In the case of a binary mixture, only the first- and 
second-degree coUisional moments have been also explic- 
itly obtained. In particular [13], 



dvmrWJrsifr, fs] = ^A*sr-(l + "rs) {2/OsPr 

Psd 



- (jrjs + jsjr) - J^P'^ril + ^rs) 
X [psPr + PrPs - (jrjs + jjr) 



+ 



- {prPs + PsPr)~ir h 



(13) 



where pr = mrrir, p = 'YlirP^ ^^'^ total mass density, 
and V = V — u is the peculiar velocity. Equations (7) 
and (8) also define the flow velocity and the partial 
temperature of species r, the latter measuring the 
mean kinetic energy of species r. Computer simulations 
[30], experiments [31] and kinetic theory calculations [32, 
33] indicate that the global granular temperature T is in 
general different from the partial temperatures Tj.. The 
mass, momentum and energy fluxes are characterized by 
the mass flux 

jr = m,ydvV/,(v), (9) 
the pressure tensor 

P = 5I / dvm,VV/,(v), (10) 



where 



P,. = / d^rmr'V^/ fr, 



(14) 



Pr = UrTr = trP^/c? is the partial pressure of species 
r, and 1 is the d x d unit tensor. It must be remarked 
that, in general beyond the linear hydrodynamic regime 
(Navier-Stokes order), the above property of the Boltz- 
mann collision operator is not sufficient to exactly solve 
the hierarchy of moment equations, due to the free- 
streaming term of the Boltzmann equation. Neverthe- 
less, there exist some particular situations (such as the 
simple shear flow problem) for which the above hierarchy 
can be recursively solved. The cooling rate C defined by 
Eq. (12) can be easily obtained from Eq. (13) as 



C - ^ XrOJrsPsri'^ + Q^rs) 



-{jrPsr + JsPrs) + 7i 



Psrji + Oirs) - 1 
dpsPr 



Jr h 



(15) 
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where 7^ = Tr/T . 



fi{V,t) reads [35] 



d 



d 



—fi - aVy — fi = Jii[/i,/i] + Jl2[/l,/2], 



dt 



(16) 



III. UNIFORM SHEAR FLOW 



We consider a binary mixture of inelastic Maxwell 
gases under USF. As mentioned in the Introduction, the 
USF state is macroscopically defined by constant densi- 
ties Hr, a spatially uniform temperature T{t) and a linear 
velocity profile u{y) — Ui{y) — U2(?/) = ay St, where a is 
the constant shear rate. Since and T are uniform, then 

= q = 0, and the transport of momentum (measured 
by the pressure tensor) is the relevant phenomenon. At 
a microscopic level, the USF is characterized by a veloc- 
ity distribution function that becomes uniform in the lo- 
cal Lagrangian frame, i.e., fr{r,v;t) ~ /r(V,t). In this 
frame, the Boltzmann equation (1) for the distribution 
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while a similar equation holds for /2. The properties of 
uniform temperature and constant densities and shear 
rate are enforced in computer simulations by applying 
the Lees-Edwards boundary conditions [37] , regardless of 
the particular interaction model considered. In the case 
of boundary conditions representing realistic plates in rel- 
ative motion, the corresponding non-equilibrium state is 
the so-called Couettc flow, where densities, temperature 
and shear rate are no longer uniform [38, 39]. 

As alluded to above, the rheological properties of the 
mixture are obtained from the pressure tensor P = Pi -|- 
P2, where the partial pressure tensors P^ (r = 1,2) are 
defined by Eq. (14). The elements of these tensors can 
be obtained by multiplying the Boltzmann equation (16) 
by m^VV and integrating over V. The result can be 
written as 



— Pi,ij + aikPi,kj 



ajkPi.ki + BiiPi^ij + Bi2P2,ij = {Aiipi + A12P2) S, 



(17) 



r 



where use has been made of Eq. (13) (with = 0). In 
Eq. (17), Oij = aSixSjy and we have introduced the coef- 
ficients 



An 



2{d + 2) 
A12 

Wii 



(l + aii)2 + 



d + 2 



M2i(l + ai2)', (18) 



^12 Pi 2 n , ^2 
7— ^— M2i(l + "12) , 
d + 2 p2 



-(1 + aii){d + 1 — an) + 



(19) 



2uj 



12 



d{d + 2)' d{d + 2) 

x^2i(l + ai2) [d + 2 - /i2i(l + ai2)] , (20) 



Bri — -'—,A\2- 
d 



(21) 



Adequate change of indices (1 o 2) provide the equa- 
tions pertaining to P2. The evolution equation for the 
temperature can be obtained from Eq. (17) and reads 



(22) 



Here, the following reduced quantities have been intro- 
duced: C* = C/i^o, a* = ajv^. P*y = Pxy/p; p = nT being 
the hydrostatic pressure. The expression for can be 
derived from Eq. (15) when one takes = 0. The result 
is 



'2 

C — "3 ^ ^ ^r*^sA^sr (l~t-^sr) 



7r 



1 + ar 



IsPrs. 

(23) 



As can be seen in Eq. (22), the temperature changes 
in time due the competition of two opposite mechanisms: 
viscous heating (shearing work) and energy dissipation 
in collisions. The reduced shear rate a* is the non- 
equilibrium relevant parameter of the USF problem since 
it measures the distance of the system from the homo- 
geneous cooling state [1]. In general, since a* does not 
depend on time, there is no steady state unless a* takes 
the specific value given by the steady-state condition 



a! P. 



s,xy 



2 



(24) 



denoting a* and Pg steady-state values of the (re- 

duced) shear rate and the pressure tensor. Beyond this 
particular case, the (reduced) shear rate and the coeffi- 
cients of restitution can be considered as independent and 
so, one can analyze the combined effect of both control 
parameters on the rheological properties of the mixture. 
This is one of the main advantages of the interaction 
model used in this paper, in contrast to previous works 
[13] : our goal is to disentangle the effects of dissipation 
from those of forcing through shear. For comparison to 
a real system or with simulation data, that are generally 
studied in the steady state where coUisional cooling and 
viscous heating equilibrate, one should however specif- 
ically work with an a-dependent (reduced) shear rate, 
given by the solution of Eq. (24) . We shall come back to 
this point in Sec. V. 

We are interested in obtaining the explicit forms of 
the scaled pressure tensors P*^^ = Pr,ii/p in the long- 
time limit. The relevant elements of these tensors are 
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P* = V 



and P*^xy ■^it'^ ^ 

As in the mono- 



id - V\P* P* 
[13, 18]. Here, p* = Pr/p = Xrjr 
component granular case [15], one can check that, af- 
ter a certain kinetic regime lasting a few collision times, 
the scaled pressure tensors P*,^ = Pr,ij/p reach well- 
defined stationary values (non-Newtonian hydrodynamic 
regime), which arc non- linear functions of the (reduced) 
shear rate a* and the coefficients of restitution. In terms 
of these scaled variables and by using matrix notation, 
Eq. (17) can be rewritten as 



1 , 2 where V is the column matrix 



V 



\ 



( P*. 

p* 

i,yy 

^l.xy 
^2,xx 
^2,yy 
\ P^.xy / 



(26) 



£V = 0, 



(25) and C is the square matrix 







d—l A * 

~a~^ii 


2a* 


^1*2 


1 A* 

d^l2 




d—1 /! + 
~3~^12 


\ 


1 A* 


Bl,+ 


\ d-l A* 







d^l2 


■^^12 


a— 1 A^ 
~d~^12 










a* 


Bli + A 












B12 


R* 1 A* 

-°21 d 21 




d-l A* 

d ^21 





^22 + 


A - ^^22 




d-l A* 
d ^22 


2a* 


1 A* 


P21 


a — 1 

~r^2i 







1 4* 

d^22 


^22 + 


\ d-l A* 
^ d ^22 





\ 







-^^21 









a* 


^2*2 + A y 



Here, A*^ = A^s/vo and B*^ = Brs/vo- In addition, 
it has been taken into account that for long times the 
temperature T(t) behaves as 

T{t) = T{0)e^''°\ (28) 

where A is also a nonlinear function of a*, a^s and the 
parameters of the mixture. The (reduced) total pressure 
tensor P*^ = Pij/p of the mixture is defined as 

P:,=PU,+Phr (29) 
Equation (25) has a nontrivial solution if 

det/: = 0. (30) 

Equation (30) is a sixth-degree polynomial equation with 
coefficients depending on a and ^ = {xi, ji, an, a22, ai2}- 
Here, fi = mi/m2 is the mass ratio. In general, this 
equation must be solved numerically. Figure 1 shows 
the real part Re[A] of the roots of Eq. (30) versus a* for 
hard disks (d = 2) in the case xi = 0.2, mi/m^ = 0.5, 
and Urs = 0.8. Obviously, exactly the same curves are 
obtained in the case mi/m2 = 2 and xi = 0.8. At a given 
value of the shear rate, the difference between the two 
largest values of i/qA gives the inverse of the relaxation 
time of the transient regime. It can be proved that this 
difference docs not vanish if xi ^ 0. 

According to Eq. (28), the largest root of the sixth- 
degree equation (30) governs the time evolution of the 
global temperature Tit) in the long-time limit. Thus, 
the upper curve in Fig. 1 gives the value of A of Eq. (28) 
for the case xi = 0.2, mi/m^ = 0.5, and a^s = 0.8. 

The stationary forms of Pi and P2 are obtained by 
solving the homogeneous equation (25) since det£ = 0. 



This equation has a nontrivial solution that can be writ- 
ten as 



v 



(31) 



where 



( Pl 

p* 

i-yy 

Plyy 
p* 

\ P2,xy I 



(32) 



6 



/Bl 



c = 



11 


' 12 12 








d 





4* 

^12 


A* 
^11 


A + Bl, 


^12 








A* 

^22 


A* 

^21 


■c'21 


A + 













a* 





A + i?i\ 












a* 


^21 


B22 + 



(33) 



A/ 



and 



/ A\^ - Bi*2 \ 
^12 



V 



4* 
^22 






(34) 



The expressions of p\ and P* ij can be obtained from Eq. 
(31). In particular, the explicit form of = 2:171 can 
be found in Appendix A. This quantity gives the ratio 
between the energy of the species 1 and the total energy 
of the mixture. 

It is important to recall that, although the scaled pres- 
sure tensors ij achieve stationary values, the binary 
mixture is not in general in a steady state since the gran- 
ular temperature changes in time. In fact, since P*y < 0, 
according to Eq. (22) the temperature T{t) grows expo- 
nentially if —2a*P*y > dC,* , namely, when the imposed 
shear rate is large enough to make the viscous heating ef- 
fect dominate over the coUisional cooling. The opposite 
occurs when d(* > —2a*P*y and so, the temperature 
decreases in time. 

To illustrate the non-Newtonian behavior [40] of the 
temperature ratio and the pressure tensor, Figs. 2 and 
3 show 71 and P*^, respectively, as functions of the (re- 
duced) shear rate a* for an equimolar mixture (xi = 0.5) 
of inelastic hard spheres (c? = 3) with mxjmi = 8. Two 
different values of the (common) coefficient of restitution 
= OL have been considered. The temperature ratio 
Ti/T measures the lack of equipartition of kinetic energy. 
As expected for driven granular mixtures [27, 30], the 
lighter particles have a smaller temperature for moder- 
ate shear rates, while the opposite happens at high shear 
rates. Figure 3 shows the dependence of the relevant el- 
ements of the total pressure tensor P*^ on a*. A signal 
of the non-Newtonian behavior is the existence of normal 
stress differences in the shear flow plane. It is also appar- 
ent that the influence of coUisional dissipation on the rhe- 
ological properties (measured through the elements P*^) 
is not quite significant, especially at high shear rates. It 
must be remarked that the trends observed for this plain 
vanilla IMM turn to be very similar to those previously 
obtained from the improved model of IMM [18] . 



IV. TRACER LIMIT (xi ^ 0) 

The results derived in the preceding section have shown 
that the time dependent solution for the second-degree 
velocity moments is given in terms of the roots of the 




FIG. 2. Shear rate dependence of the temperature ratio Ti/T 
in three dimensions (d = 3), for an equimolar mixture (xi = 
0.5) with m\/m2 — 8. Two values of the (common) coefficient 
of restitution a^s = a have been considered: a = 1 (solid line) 
and a = 0.8 (dashed-dotted line). 



sixth-degree polynomial equation (30). For long times, 
the dominant behavior is described by the two real roots, 
Ai and A2. In particular, the energy ratio Ei/E = xi'yi 
can be written as 



Ph ^ Apl{\2 ,a,i) 
E 



-(A2-Ai)iyot 



(35) 



where A and B are constants depending on the initial 
conditions and the explicit expression of p* is given by 
Eq. (Al). After a relaxation time of the order of |(A2 — 
Ai)t'o|^^, the energy ratio Ei/E reaches a steady state 
value pt(Aniax,a,l) where A^ax = max(A2, Ai). As long 
as xi 7^ 0, one has Ai 7^ A2 for any value of the shear 
rate and ^. 

Let us assume now that the mole fraction of one of the 
species (say for instance, species 1) becomes negligible. 
In the tracer limit {xi — > 0), the sixth-degree equation 
(30) for A factorizes into two cubic equations (see ap- 
pendix B, where the quantities used below are defined), 
with the following largest roots: 



X^^^^2F{ayAg^)Ag^- 



,(0) 

^22 



l(0) 
^22 



where 



(1 + 0^22)^ 
d + 2 



Fix) = - sinh^ 
^ ' 3 



F{Z) 



*22 



2d 



^ cosh ^ ( 1 + ^.T^ 



(36) 



(37) 
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and 




0.0 0.2 0.4 0.6 0.8 1.0 



a* 

FIG. 3. Shear rate dependence of the (reduced) elements of 
the pressure tensor P^^ (a), Pyy (b) and P^y (c). Two values 
of the (common) coefRcient of restitution ars = a have been 
considered: a = 1 (solid lines) and a — 0.8 (dashed lines). 



Xi{a*,Xi)^xf\a*) + \['\a*)xu (42) 

where A2''^ and A^^"* are given by Eqs. (36) and (39), re- 
spectively. The expressions of A2^^ and A^^^ can be ob- 
tained from the general sixth-degree polynomial equation 
(30); they are given in Appendix C. It can be checked 

that if A = A^°' in Eq. (40), then [according to Eqs. (CI) 
and (C2)] D{xf^) ^ 0, Ao(A^°^) ^ 0, and so the energy 
ratio Ei/E vanishes when xi — > 0, as may have been 
expected. However, if A = xf^ in Eq. (40), Ao(A^"^) = 
which implies that Ei/E ^ 0. Therefore, by taking the 
tracer limit in Eq. (40), one gets 

r 1 *a r * ^ * ^ D{xf(a*),a*) 

hm — A2(a ,a;i),a ,xi = -, (43) 

^1^0 2^1 Ao(A^"^(a*),a*) 



and 



a* _ 2{d+2) 
If" (l + a22y 



(38) 



lim p*(Ai(a*,.Ti),a*,a:i) 

!-0 



D{X^^\a*),a*) 



Aoi(a*) + Ai(A(°Ha*),a*) 
^ 0, (44) 



The above root A2'''' rules the dynamics of the host fluid 
(excess component) while the evolution of the impurity 
is governed by the root 



(0) 



2F{a*/Af^)Af^ 



d + 2 
2 



(l + ai2)"F 



/ a (1-Ha22)^ 

Wl, (1 + «12)2 
Ai21 



where 



Aoi(a*) 



aAo(A,a*) 
dX 



X['\a*). (45) 



Ai2i(l + ai2) 1-— (1 + 012) ■ (39) fact that Ao(Ai' 



Equation (43), where A2 is the argument oipi, is relevant 
for the case Aj^'' > A^^'' while conversely, Eq. (44) applies 
when A2"'' < X^i \ In Eq. (44), use has been made of the 



0. 



As seen above, the largest of all roots, Amax, is the rele- 
vant one to obtain the asymptotic energy ratio Ei/E. Wc 
now show that the behavior of the system is qualitatively 



very different depending on Amax = A^^*'' or Amax = A2 . 
For xi 7^ 0, the expression of pi (that corresponds to the 
long-time value reached by Ei/E) is given in Appendix 
A by Eq. (Al). If xi — ;> 0, the energy ratio pi becomes 



(0) 



P*ii\a*,$,) « xi 



D{X,a*) 



Ao(A,a*)+Ai(A,a*)a;i' 



(40) 



where the dependence on ^, an, a22, and ai2 is implic- 
itly assumed on the right-hand side. The general expres- 
sions of D, Aq, and Ai can be found in Appendix C. 

Equation (40) holds for both A = A2 and A = Ai . These 
two possibilities turn out to differ in that A^*^' (the value 
of Al at xi = 0) is a root of Aq. It is then important 

to keep track of the finite xi correction to A^*^' that is 
present in Ai. To first order in xi, we have 



A2(a*,xi) 



Xf\a*) + Xi'\a*)xi 



(41) 



In conclusion, if A,, 



= X^\ the right-hand side of 
Eq. (43) is the temperature ratio Ti/T and so, the en- 
ergy ratio Ei/E = 0. On the other hand, if Amax = xf\ 
the temperature ratio diverges to infinity and the en- 
ergy ratio becomes finite. The change from one behavior 
to the other is akin to a non-equilibrium transition be- 
tween disordered {Ei/E — 0, finite temperature ratio) 
and ordered {Ei/E ^ 0, diverging temperature ratio) 
phases [23]. The task now boils down to identifying the 

regions of parameter space where a''''' = Aj''^ from which 
the domains of existence of the ordered and disordered 
phases can be obtained. It is important to note here 
that although our procedure leads to tracer limit quan- 
tities that a priori depend on the intruder-intruder coef- 
ficient of restitution an, such a parameter turns out to 
disappear from the final expressions, which is intuitively 
expected. We could not analytically show this property 
(the computer package of symbolic calculations (MATH- 
EMATICA) used was not able to provide a manageable 
expression), that is nevertheless a systematic numerical 
observation. We will henceforth drop an from the rele- 
vant parameters in the tracer limit. 
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A. Absence of shear rate (homogeneous cooUng 
state) 

Before analyzing the physical consequences of the 
above niathematical treatment in the shear flow prob- 
lem (a 7^ 0), it is instructive to consider the particular 
case of vanishing shear rates. This state is referred to as 
the homogeneous cooling state (HCS). Such a situation 
has been analyzed in detail in a recent work [41]. Here, 
for the sake of completeness, we summarize the most rel- 
evant results in the HCS. 



When a* = 0, the roots X^2^ and A^"'' simply reduce to 



(0) 



AO) 



ah 



2d 



AO) 



; ^2l(l + Ctu) 



M21 



(1 + ai2) 



(46) 
(47) 



It is apparent that, even in the HCS, there are two dif- 
ferent regimes of behavior depending if A2'^'' is smaller 



or larger than 
critical mass ratios 



Equating x'^^ 



and a1°^ 



leads to two 



ai2 



l+ai 



and 



(+) _ 



l+ai2 
(48) 



with /i[j(]s < ^[j+'g. When /x[jjg < ^ < ^[j+^g, we have 



Aj'^'' > A^'^'' and the temperature ratio Ti/T remains finite 
(disordered phase). It is given by [41] 



v(0) 



2Atl2M2l(l + Ckl2)^ 



T 4/i2i(l + ai2) [1- ^(l + ai2)] - 1 



(49) 



The expressions of /i^j^g and Ti /T coincide with the ones 
previously derived by Ben-Naim and Krapivsky [28] in 
their analysis on the velocity statistics of an impurity 
immersed in a uniform granular fluid. When /.j > ^^q^ 

or fi < fi\[(lg, Xf^ < Xf\ so that Ti/T ^ oo but the 
energy ratio Ei/E is finite. This is a new result together 

with the identification of the bound /^hcs- '^^^ explicit 
expression of Ei/E is [41] 



El 

E 



- 1 + A^J^2l{l + ai2) [1 - ^(1 + ai2)] 



^22 



l + 2M2i(l-a22) 



(50) 



Three remarks are in order here: 

• while the upper bound /^[j^g is well defined for all 
values of ai2 and 0:22, the lower one is meaning- 
ful (i.e. positive) only when ai2 > \/{l + o^2)/^- 
Such a constraint cannot be met when ai2 = a22 
nor when ai2 < •\/2/2, and requires "asymmetric" 
coefficients of restitution (the above inequality im- 
plies ai2 > 0:22, i.e. more dissipative inter host gas 
collisions than cross intruder-gas encounters). 



Ordered phase 
^ light tracer 



= 0.99 
= 0.5 

= 



Disordered phase 



Ordered phase 
lieavy tracer 



,(+) 



. 1000 

,.(+) 



FIG. 4. Phase boundaries discriminating, in the tracer limit 
xi — >■ and for hard splieres {d — 3), between ordered and 
disordered phases. In the ordered phases, the energy ratio 
El /E (impurity over total energy) is finite -hence a diverging 
temperature ratio Ti/T- whereas in the disordered phase, the 
order parameter Ei/E vanishes with a corresponding finite 
temperature ratio. The vertical arrows in the lower panel 



indicate the threshold f^^^Qs as given by Eq. (48). The values 
of ^[[^' — \/2—l [here, common to all three parameter sets, see 
Eq. (51)] and of a*'"^' [given by Eq. (53)] are also indicated. 



• a similar extreme breakdown of the energy equipar- 
tition has been found [42] for inelastic hard spheres 
since, in the ordered phase [Ti/T — > 00). the ra- 
tio of the mean square velocities for the impurity 
and gas particles m2Ti / miT becomes finite for an 
extremely heavy impurity particle (toi/to2 00). 
Our Maxwell approach hence captures robust ef- 
fects, but at the same time exaggerates the trend 
of more refined models. 

• in the present unforced situation, the results are 
independent of the space dimension d. 



B. Nonzero shear rate 

We now wish to assess the influence of the shear rate 
on the transition observed in the absence of shear (HCS). 
When a* = 0, we have seen that an ordered phase ex- 
ists for heavy impurities, while asymmetric dissipation 
may open a window for a light impurity ordered phase, 
provided ai2 > 1/(1 -I- ah)/'^- On the other hand, it 
is known for elastic systems that an ordered phase sets 
in for jj, < \/2 — 1 [43] (which can be seen as a light 
impurity condition), provided the shear rate a* is larger 
than a certain critical value a*(/x). To see how these 
two limiting cases are connected, we first show in Fig. 
4 how the shear rate affects the HCS scenario, and how 
coUisional dissipation modifies the results obtained for 
ordinary gases (a^s = 1) [23]. Since we have chosen sym- 
metric dissipation parameters {ai2 = 0:22) in Fig. 4, the 
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analytical results. In particular, the threshold ^[j^ may 



FIG. 5. Same as Fig. 4, for liard disks (d = 2). Here, Q12 = 
a-21 = 0.9. 



light tracer ordered phase is precluded on the a* = axis. 
It can be seen that this phase exists, in the shear rate 
versus mass ratio plane, provided a* > a* and /i < 
which defines a domain that is rather insensitive to the 
value of the coefficient of restitution. On the other hand, 
the heavy tracer ordered phase is much more sensitive to 
coUisional dissipation. At variance with the light tracer 
phase, it has an enhanced domain of existence for more 
dissipative systems, and disappears in the elastic limit, 
as suggested by Fig. 4 (see the dashed line corresponding 
to ai2 = a22 = 0.99, squeezed in the lower right corner 
of the graph) . We also note that the heavy impurity or- 
dered phase is destroyed by a sufficiently vigorous shear 
rate, see below. 

The characteristic features seen in Figure 4 -with two 
distinct ordering pockets- may be rationalized by further 
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= 0.9 ; = 0.9 
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Disordered phase 



Ordered phase 
heavy tracer 



t 



^^HCS 



FIG. 6. Tracer limit phase diagram, in three dimensions. The 
values of ^^^'3 and /^^'s predicted by Eq. (48) are shown by 
the arrows. For Q12 = 0.9 as chosen here, the constraint ai2 > 
•\/(l + 022)72 for the existence of a light impurity ordering at 
vanishing shear, reads 022 < 0.787. The dashed line therefore 
does not reveal any light impurity ordered phase at zero shear 
rate (this was also the case in Fig. 4). 



be obtained by the condition \\ 
Following that route, we find 



(0) 



1 + ai2 
1 + a22 



aI*^^ when a* — > 00. 



(51) 



that does not depend on space dimension, as the thresh- 
olds pertaining to the HCS. As expected, we recover the 
threshold value = \/2 — 1 for clastic systems. Like- 
wise, the upper shear rate a*'^^ beyond which the heavy 
impurity ordered phase disappears may be derived from 
enforcing the limit yu 



0, then 



00 in the equality Aj^' = A^''^ 



Since ^21 

implies [see Eq. (36) 
F(a) 

Equation (52) yields 



so that A 



(0) 



d + 2l 



"22 



2d 1 



0. This 



(52) 



!*(+) 



\ + d — OL22 \ 1 



2(d+2)' 



where use has been made of the relation 



X 

~d 



(53) 



(54) 



where F{x) is given by Eq. (37). The identity (54) allows 
for a convenient expression of a [and hence a* through 
Eq. (38)], once A [and hence F(a)] is known. Since 
the value of a**^"*"^ has been obtained from the condition 
Aj"'' = (constant temperature in the ordered phase), 
the expression (53) also gives the a-dependence of the 
shear rate in the steady USE state [18]. 

The quantities A^th^ ^'^'^ a*^"*"-* are indicated by arrows 
in Fig. 4. Consistent with the disappearance of the heavy 
impurity ordered phase for elastic systems, is the vanish- 
ing of a*^"'"-' when 022 — >■ 1- It can also be noted that 
is only a host property, and does not depend on 
a\2- Finally, Fig. 5 shows that the behavior in two and 
three dimensions are similar. 

We now turn to asymmetric coUisional dissipation 
cases, so that the light impurity ordered phase may exist 
even for vanishing shear rates. Such a scenario is illus- 
trated in Fig. 6, that corroborates the analytical predic- 
tions. The boundary of the light impurity ordered phase 
in a shear-mass ratio phase diagram is non trivial, and 
indicates the existence of an interval of /i values, below 
/^HCS' '^ith 8- re-entrance feature. Indeed, starting from 
the ordered phase at a* = 0, and increasing a* at fixed 
/i, one first meets a transition from order to disorder, 
followed by a subseque nt ordering t ransition. Similarly, 
and again for a\2 > \/(l~-i-"a^^)/2, the following series 
order — > disorder — > order occurs when fj, is increased at 
fixed reduced shear rate, provided a* < a*^~^\ 
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FIG. 7. Order parameter E-i/E as a function of the reduced 
shear rate a*, for different values of q = ai2 = 022- The main 
graph is for ^ = 0.1 (light impurity). In the inset, where the 
results for /i = 20 are displayed with the same convention for 
the different curves as in the main graph, the value of a = 0.9 
is not shown, since with such a parameter, /i^'g > A^i hence 
no ordering is possible. In other words, the counterpart of the 
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FIG. 8. Order parameter for d = 3, Q12 = 0.9 and 022 = 0.55 
(see the dotted line in Fig. 6). The value jjL = 10~^ leads to a 
disordered phase at zero shear rate (thick continuous curve), 
while the smaller values of reported here are associated with 
small shear ordering, with the possibility of an intermediate 
disordered phase (see text). 



To substantiate the phase diagram reported above, we 
show in Fig. 7 the shear rate dependence of the order pa- 
rameter Ei/E, for different values of the (common) coef- 
ficient of restitution and mass ratios. The light impurity 
ordering is seen to be enhanced by increasing the shear 
rate, while the reverse behavior is in general observed for 
heavy impurities (see the inset). The critical thresholds 
observed in Fig. 7 are fully compatible with those ap- 
pearing in Fig. 4. Focussing next on the light impurity 
ordered phase, we report the order parameter variation in 
cases of asymmetric coUisional dissipation. To this end, 
we return to the set ai2 = 0.9 and 022 = 0.55 addressed 
in Fig. 7. For such quantities, one has /^ncg — 0.051. 



FIG. 9. Shear rate dependence of Ei/E for hard spheres 
(d = 3) in the case fi = 0.2 and 0:22 = 0:12 = 0.9. Three 
different values of the mole fraction xi have been considered: 
XI — 0.1 (circles), xi — 0.01 (squares), and a;i = (solid 
line). 



As can be seen in Fig. 8, when ^ < /x^Jg, the order 
parameter is non vanishing at small shear rates, while 
when /Xygg < /.j < fi[-^ , ordering sets in only beyond a 

critical shear rate (note that ju[^'' ~ 0.733, so that the 
light impurity ordered phase does exist in some shear do- 
main, for all the mass ratios used in Fig. 8). The figure 
clearly shows the re-entrance of order alluded to above 
(see the curves for ^ = 10~^ and /i = 3xl0~^, where an 
intermediate interval of a* values leads to a disordered 
phase with Ei/E = 0. For ^ = lO"-^ (dashed-dotted 
line), all shear rates lead to phase ordering {Ei/E 7^ 0), 
but there is a fingerprint of the reentrant behavior in the 
non monotonicity of the order parameter with a*. 

For completeness, we now show how the tracer limit 
is approached. For /x < and ai2 = 0:22, order- 

ing occurs when the system is driven sufficiently far from 
equilibrium (a* > a*). To illustrate how the abrupt tran- 
sition observed in the tracer limit is blurred by finite con- 
centration, Fig. 9 shows Ei/E versus the (reduced) shear 
rate a* for d = 3, mi/m2 = 0.2, an = 0.21 = ai2 = 0.9 
and X\ = 10~^, 10~^ and 0. It is apparent that the curves 
tend to collapse to the exact tracer limit result as the 
mole fraction x\ vanishes. This is indicative of the con- 
sistency of the analytical results derived at X\ ~ 0. More- 
over, since the impurity particle is sufficiently lighter 

than the particles of the gas (/x = 0.2 < [i'^ ~ 0.414), 
the energy ratio E\ / E is non vanishing in the tracer limit 
if a* > a* ~ 7.557 in the present case (see the figure). It 
can also be noted that at finite xi, the energy ratio for 
small shear rates is of the order of xi, since the temper- 
ature ratio is quite close to unity in that limit (see e.g. 
Fig. 2). 

Although we have focused our attention on the energy 
ratio, it is clear that similar features can be analyzed 
when considering other quantities. An interesting candi- 
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FIG. 10. Shear rate dependence of the intrinsic viscosity [77*] 
for hard spheres (d = 3) in the case mi/m2 — 0.1. Three 
different values of the (common) coefficient of restitution a = 
ai2 = a22 have been considered: a = 0.5 (soHd hne), a — 0.7 
(dashed line), and a — 1 (dotted line). 



date is the non-linear shear viscosity 77* defined as 



(55) 



where the shear stress P*y = P^^xy + P^^y The rheo- 
logical function 77* characterizes the nonUnear response 
of the system to the action of strong shearing. In terms 
of p* and for finite values of xi, the expressions of Pi ^y 
and P2,xy ^'^^ given by Eqs. (A6) and (A7), respectively. 

As expected, in the tracer limit [xi — > 0) and in the 
disordered phase, the total shear viscosity 77* of the mix- 
ture coincides with that of the solvent gas 77* — > rj*. where 
77* is given by Eq. (C18). However, in the ordered phase, 
there is a finite contribution to the total shear viscosity 
coming from the tracer particles (see Eq. (C19)). To il- 
lustrate it, Fig. 10 shows the shear rate dependence of the 
intrinsic shear viscosity [44] [77*] = \mix^^o{vi* — r]s)/vt 
for the mass ratio mi/m^ = 0.1 (hght impurity) and dif- 
ferent values of the (common) coefficient of restitution. 
We observe that the intrinsic viscosity [77*] is clearly dif- 
ferent from zero for shear rates larger than its correspond- 
ing critical value. On the other hand, its magnitude is 
smaller than the one obtained for the order parameter 
Ei/E (see Fig. 7) 



V. DISCUSSION AND CONCLUSION 

In this paper we have analyzed the dynamics of an im- 
purity immersed in a granular gas subject to USF. The 
study has been performed in two successive steps. First, 
the pressure tensor of a granular binary mixture of in- 
elastic Maxwell gases under USF has been obtained from 
an exact solution of the Boltzmann equation. This solu- 
tion applies for arbitrary values of the shear rate a and 
the parameters of the mixture, namely, the mole fraction 



xi = ni/n, the mass ratio fi = mi/m2 and the coeffi- 
cients of restitution an, a22 and ai2. Then, the tracer 
limit {xi 0) of the above solution has been carefully 
considered, showing that the relative contribution of the 
tracer species to the total properties of the mixture does 
not necessarily vanish as xi — > 0. This surprising result 
extends to inelastic gases some results derived some time 
ago for ordinary gases [23]. The above phenomenon can 
be seen as a non-equilibrium phase transition, where the 
relative contribution of the impurity to the total kinetic 
energy Ei/E plays the role of an order parameter [45] . 

The transition problem addressed here has been ana- 
lyzed in the framework of the Boltzmann equation with 
Maxwell kernel. The key advantage of inelastic Maxwell 
models, in comparison with the more realistic inelastic 
hard sphere model, is that the collisional moments of the 
Boltzmann collision operators Jrs[fnfs\ can be exactly 
evaluated in terms of the velocity moments of fr and /s, 
without the explicit knowledge of these velocity distribu- 
tion functions. Here, we have explicitly determined the 
collisional moments associated with the second-degree ve- 
locity moments to get the pressure tensor of the mixture. 
In addition, the collision rates Urs appearing in the op- 
erators Jrs[fr, Is] have been chosen to be time indepen- 
dent so that the interaction model allows to disentangle 
the effects of collisional dissipation (accounted for by the 
coefficients of restitution) from those of boundary condi- 
tions (embodied in the reduced shear rate a* defined as 
a* ~ a/uQ, Uq being a characteristic collision frequency). 
Consequently, within this model, collisional dissipation 
and viscous heating generally do not compensate, so that 
the granular temperature increases (decreases) with time 
if viscous heating is larger (smaller) than collisional cool- 
ing. 

In our system, the temperature T2 of the gas particles 
(T2 ~ r in the disordered phase when xi — >■ 0) changes in 
time due to two competing effects: the viscous heating 
term {—aPxy) and the inelastic cooling term (C?2). In 

fact, for long times, T2{t) ~ exjp[\^2^ voA where Aj*^' = 
-C* - {2a*/d)P*y (here, C* = C/^^o and = P^y/p). 
Since P*y < 0, the cooling rate can be interpreted as 
the "thermostat" parameter needed to get a stationary 
value for the temperature T2 of the gas particles. At a 
given value of the shear rate, the cooling rate increases 
with dissipation and so, T2 decreases as Urs decreases. 
The tracer particles are also subject to two antagonistic 
mechanisms. On the one hand, Ti — >■ 00 due to viscous 
heating and on the other hand, collisions with the gas 
particles tend to "thermalize" Ti to T2. Both effects are 
accounted for by the root A^"'' that governs the behavior 
of Ti for long times, i.e., Ti(t) ~ exp[A^°'t'o<]. When 



> X'^\ the temperature ratio T1/T2 grows without 



bounds. The parameter ranges where such a requirement 
is met define the ordered "pockets" of the phase diagram, 
and where the energy ratio Ei / E -explicitly worked out 
here- reaches a finite value. We have found that two 
different families of ordered phase can be encountered 



• a light impurity phase, provided that the mass ratio 
/i does not exceed the threshold ^[^^ given by Eq. 
(51). Such a phase always exists for shear rates 
larger than a certain critical value, but can also be 
observed at vanishing shear in cases of asymmetric 
collisional dissipation, whenever gas-gas collisions 
are sufRcicntly more dissipative that intruder-gas 
collisions (ai2 > y^l + Q!|^)/2 ) . 

• a heavy impurity phase, that -unlike the light im- 
purity phase- cannot accommodate large shear, 
and requires a* < a*^'^\ where the threshold 0*^+^ 
is given by Eq. (53). 

The fact that a* and a^s are independent parameters 
[unless a* takes the specific value a* given by the steady 
state condition (24)] allows one to carry out a clean an- 
alytical study of the combined effect of both control pa- 
rameters on the properties of the impurity particle. It is 
however important to bring to the fore the precise cou- 
pling between shear and collisional dissipation that the 
steady state condition Eq. (24) implies. The answer de- 
pends on the phase considered, ordered or disordered and 
the energy balance embodied in Eq (24) can be expedi- 
ently expressed as max(A^°\ Aj^'') = 0. In other words, 
a* follows from enforcing x'^^ = in the disordered state, 

and likewise X^f^ = in the ordered regime. The solu- 
tion g of the first equation has already been displayed 

in Eq. (53) : 02 .s a*^^'. On the other hand, the solu- 
tion at to aJ°^ = can be obtained along similar lines 
and reads 

• ^3/2 / l-^/?i2 _ d + 2 f 2 \] 

where we have introduced the notation /3i2 = /i2i(l + 
ai2). Finally, the steady state condition implies a* ~ 
(resp. a* = g) in the disordered (resp ordered) 
case. The corresponding line in the shear versus mass 
ratio plane is shown in Fig. 11, for a given set of dis- 
sipation parameters that corresponds to one of the sit- 
uations analyzed in Fig. 6. It appears that remaining 
on the steady state line shown by the thick curve, one 
spans the three possible regimes: light tracer ordering 
at small /x, disordered phase at intermediate values, and 
heavy tracer order at larger fi. We conclude here that 
the scenario uncovered in our analysis is not an artifact 
of having decoupled shear from dissipation, and we em- 
phasize that from a practical point of view, studying the 
transient regime (before the steady state occurs) anyway 
offers the possibility to enforce the above decoupling. 

In addition, we would like to stress here that, in spite 
of the approximate nature of our plain vanilla Maxwell 
model, the results obtained for binary mixtures under 
steady USF [13] compare quite well with Monte Carlo 
simulations of inelastic hard spheres [19]. This can be 
seen in Fig. 12 for the (reduced) elements of the pressure 
tensor P* in the steady shear flow state defined by the 




FIG. 11. Same plot as Fig. 6, where only the curves per- 
taining to the values ai2 ~ 0.9 and Q22 = 0.55 have been 
retained. The ordered pockets are shown by the hatched ar- 
eas. The steady state condition (24) of equal viscous heating 
and coUisional dissipation, implies a coupling between a* and 
/i, that corresponds to the thick continuous line. The bell 
shaped dashed curve shows the locus of points where Aj"' = 0. 
It therefore coincides with the thick continuous hne in the or- 
dered phase. However, it does not yield the steady state con- 
dition in the disordered case, where the reduced shear rate is 
given by 02^3 = a*^'^\ that does not depend on /i (hence the 
horizontal sector in the thick continuous hne). 



condition (24). Therefore, we expect that the transition 
found in this paper is not artefactual and can be detected 
in the case of hard spheres interaction. 

It is quite natural -when analyzing the dynamics of an 
impurity immersed in a background of mechanically dif- 
ferent particles- to invoke two assumptions. First, that 
the state of the solvent (excess component 2) is not af- 
fected by the presence of the tracer (solute) particles 1. 
Second, that the effect on the state of the solute due to 
collisions among the tracer particles themselves can be 
neglected. We have seen that the second expectation is 
correct (the coefficient of restitution an is immaterial in 
the tracer limit, a property that is not obvious from the 
cumbersome analytical formulas reported here), but the 
first expectation is invalidated in the regions where the 
ordered phase sets in. Consequently, the seemingly nat- 
ural "Boltzmann-Lorentz" point of view -with the one- 
particle velocity distribution function /2 of the granular 
gas obeying a (closed) nonlinear Boltzmann kinetic equa- 
tion while the one-particle velocity distribution function 
/i of the impurity particle obeys a linear Boltzmann- 
Lorentz kinetic equation- breaks down. Our results show 
that collisions of type 2-1 affect /2, despite being much 
less frequent than collisions of type 2-2. We conclude here 
that, rather unexpectedly, the tracer problem is as com- 
plex as the general case of a binary mixture at arbitrary 
mole fractions. 
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FIG. 12. (color online) Plot of the reduced elements of the 
pressure tensor as functions of the (common) coefficient of 
restitution a — an — Q12 — Q22 for a three dimensional 
system (d = 3) in the steady USF. The predictions of the 
inelastic Maxwell model (solid line, present work) are tested 
against Monte Carlo simulation data (symbols, taken from 
Ref. [19]). The parameters are xi — 0.5 (equimolar mixture) 
and fi = 2. It should be noted that for a meaningful com- 
parison, a* and a need to be coupled (see Eq. (24)), since 
inelastic hard spheres enjoy a steady state for a precise value 
of the reduced shear rate, that depends on the coefficient of 
restitution a. 
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Appendix A: Energy ratio p* and shear stress P^y 

The expression of the energy ratio pi can be obtained 
from Eq. (31). It can be written as 
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S = d{Al,-Al2-Bl,+Bl.,-\) 

X [A^ + [Bl, + 522)A + 5*^522 - Bl^Bli] . 

(A5) 

Once the energy ratio is known, the remaining relevant 
elements of the pressure tensor can be easily obtained 
from Eqs. (31)-(34). In particular, the shear stress P*y 
is given by P*y = Pl -^y + Pl^y where 

d [AI2 - 5i*2 - (5i*i +X-Al,+ AI2 - Bl2)pl] 
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2,xy 
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Appendix B: Derivation of Aj'^' and 



When Xi — >■ 0, the sixth-degree equation (30) for the 
rates A factorizes into two cubic equations given by 
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where aI^J and 5rs'' denote the zcroth-order contribu- 
tions to the expansion of A*^ = Ars/i^o and B*^ = 
Brs/^o, respectively, in powers of xi. They are given 
by 
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Equation (Bl) is associated with the time evolution of 
the excess component. Its largest root is given by Eq. 
(36). On the other hand, Eq. (B2) gives the transient 
behavior of the impurity. Its largest root is given by Eq. 
(39). 



Appendix C: Some explicit expressions in the tracer 

L = d(5i*2-At2) [A^ + {Bl, + B*2)\ + Bl,B*2 - Bl^B*,] ' , li™* 

(A3) 
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In this Appendix, we provide some of the expressions used along the text in the tracer hmit. First, the quantities 
Z?, Aq and Ai appearing in Eq. (40) are given by 
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In these equations, A^^J and S^r'' arc given by Eqs. (B3)-(B6) and 
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The expressions of the coefficients xf'\a) can be written as 
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Here, we have introduced the quantities 
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where 
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With respect to the shear stress P* , in the disordered 



phase (A2°^ > A^'"'' and so, pi = 0) one has Pj* ^.^ = 
while Pi^xy is given by (see Eq. (A7)) 



k(O) 



p* 

^2,xv - 



A. 



(0) 



(B 



(0) 

22 



Af) 



(C15) 



p* _ 



4(0) _ n(0) _ 
^21 -°21 



B(°^ + Ar-A 



(0) 



l(0) 



(0) 
21 



B. 



(C17) 



where pi is given by Eq. (44). Consequently, according 
to Eq. (55), the non-linear shear viscosity jy* in the dis- 
ordered phase is 
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while in the ordered phase the result is 
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In the ordered phase (A^^"* > A2''' and so, pi ~ finite), 
Eqs. (A6) and (A7) yield, respectively. 
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